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This work reviews recent advances in the analytical treatment of the 
continuum spectrum of correlated few-body non-relativistic Coulomb 
systems. The exactly solvable two-body problem serves as an introduc- 
tion to the non-separable three-particle system. For the latter case we 
discuss the existence of an approximate separability of the long and the 
short-range dynamics which is exposed in an appropriately chosen curvi- 
linear coordinates. The three-body wave functions of the long-ranged 
part of the Hamiltonian are derived and methods are presented to ac- 
count approximately for the short-ranged dynamics. Furthermore, we 
present a generalization of the methods employed for the derivation of 
the three-body wave functions to the scattering states of iV charged 
particles. To deal with thermodynamic properties of finite systems we 
develop and discuss a recent Green function methodology designed for 
the non-perturbative regime. In addition, we give a brief account on 
how thermodynamic properties and critical phenomena can be exposed 
in finite interacting systems. 



1. Introduction 

The theoretical treatment of many-body systems is of a fundamental importance 
for a variety of branches in physics. Of particular relevance to this work are highly 
excited systems that consist of a finite number of charged particles. Such systems 
are encountered in processes involving simultaneous excitations of few particles, e.g. 
as is the case in a multiple ionization reaction. From a theoretical point of view it 
is fortunate that the scattering states of two particles interacting via the Coulomb 
potential are known exactly, for such states can serve as a benchmark for approxi- 
mate methods that are developed to deal with more complex systems. For example, 
the perturbative (Born) treatment of the two-particle Coulomb scattering and the 



comparison with the exact results gives a first hint on the problems connected with 
treatment of Coulomb scattering: The first order term of the Born perturbation 
series delivers already the correct results for the scattering cross section. Higher 
order terms are however divergent. This problem is mainly traced back to the infi- 
nite range of the Coulomb interaction that prohibits free motion in the asymptotic 
regime. 

In the theoretical treatment of more than two particle systems a further funda- 
mental difficulty arises which stems from the inherent non-separability of many-body 
interacting systems. This problem is of a general nature and appears basically for 
almost all forms of the inter-particle interactions. 

This work provides an overview on recent theoretical efforts to tackle analytically 
the problems connected with the infinite-range tail of Coulomb potentials and the 
non-seprabale aspects of many-body systems. As will be shown below, generally 
this can be done only in an approximate way. To deal with the long-range Coulomb 
interaction one dismisses the use of standard many-body approaches and attempts 
at a direct (approximate) solution of the many-body non-relativistic Schrodinger 
equation. To expose the general feature of Coulomb potential scattering and to 
introduce the basic methods used to solve the Schrodinger equation we start with 
a brief review of the two-particle Rutherford scattering. The three-body Coulomb 
scattering problem has received recently much of attention |2Hl3|. In this work we 



focus only on one aspect of this problem, namely the existence of an approximate 
separability that allows to derive three-body wave functions valid in certain region 
of the Hilbert space. For the general case of a system consisting of N interacting 
charged particles we derive correlated scattering states that are to a first order in 
the distance exact at larger inter-particle separation. 

The wave function approach is of a limited value when it comes to the study 
of thermodynamic properties and critical phenomena in finite systems. In this case 
detailed information on the density of states is needed. Such information is encom- 
passed in the many-body Green function. Therefore we devote a section of this work 
to a general scheme for the derivation of the many-body Green operator and briefly 
review a possible method to extract thermodynamic information on finite systems 
starting from the Green function. 

Unless otherwise stated we employ atomic units throughout and neglect rela- 
tivists corrections. 



2. Two charged particle scattering 

To introduce the general frame work and the basic notation let us consider the non- 
relativistic scattering states of two charged particles with charges z\ and z^. The 
Schrodinger equation describing the motion in the two-particle relative coordinate 
r is 



2/i r 



*k(r) = 0. (1) 



Here k is the momentum conjugate to r and E = k 2 /2fi is the energy of the relative 
motion, /x is the reduced mass of the two particles. To decouple kinematics from 
dynamics we make the ansatz: 



# k (r) = e ikr * k (r). 



(2) 



The distortion factor ^ in Eq.(0) is solely due to the presence of the potential. The 
asymptotic properties of (|l|) follows upon substitution of (fj) in (|l]). Then, terms 
that fall off faster than the Coulomb potential can be neglected which leads to the 
equation 
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*k(r) = 0. (3) 



This equation can be solved by the ansatz \l/ = exp(i0). Upon insertion in Eq.( 
this ansatz yields 
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W)-->-—r- Ina(r^k-r). (4) 



The factor z\Zi\ijk is called the Sommerfeld parameter and characterizes the strength 
of the interaction. The integration constant a has a dimension of a reciprocal length 
and the value a = k. The important point here is that the natural coordinate that 
appears in the treatment of Coulomb scattering is the so-called parabolic coordi- 
nate £ := r =F k • r where the + or — corresponds respectively to incoming or 
outgoing-wave boundary conditions. 



3. The three-particle coulomb continuum states 

In contrast to the two-body problem, an exact derivation of the three-body quantum 
states is not possible. Nonetheless, under certain (asymptotic) assumptions analyt- 
ical solutions can be obtained that contain some general features of the two-body 
scattering, such as the characteristic asymptotic phases. As in the preceding section 
the center-of-mass motion of a three-body system can be factored out. The internal 
motion of the three charged particles with masses raj and charges Zj ; i £ 1, 2, 3 can be 
described by one set of the three Jacobi coordinates (r^, R&); i, j, k G {1, 2, 3}; ey* ^ 
0; j > i. Here ry is the relative internal separation of the pair ij and R& is the po- 
sition of the third particle (k) with respect to the center of mass of the pair ij. The 
three sets of Jacobi coordinates are connected with each other via the transformation 

SHU:) - d (r:H(r:) < 5 > 



where 
Di = 
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D 1= ( , -/Wm 2 1 \ (6) 
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The reduced masses are defined as % = mivajjimi + nij) ; i,j € {1,2,3}; j > z. 

Accordingly, the momenta conjugate to (ry,Rjt) are defined as (k^K^). These 
momenta are related to each other by 

k23 \ t\ t ( ^12 \ t\ t ( ki3 



where D/ and D 2 * are transposed matrices of Dx and D 2 , respectively. The scalar 
product (ry, Rfc) • ( j is invariant for all three sets of Jacobi coordinates. The 
kinetic energy operator Hq is then diagonal and reads 

#o = ~^-A ri] - -^A Rk V (r y , R fc ) , (8) 
where = m^{rrii + m 3 ) / (mi + m 2 + m.3). The eigenenergy of (|8]) is then given as 

E = ^ + ^- V (ry.H*). (9) 
The time-independent Schrodinger equation of the system reads 
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H0 + -E&-E 
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(r w ,R m |* kfc!jKm > = 0. (10) 



Here we defined the product charges z y - = ZiZf, j > i G {1,2,3}. The relative 
coordinates r y - occurring in the Coulomb potentials have to be expressed in terms 
of the appropriately chosen Jacobi-coordinate set (r^,R m ). 

Asymptotic scattering solutions of (|T0D for large interparticle distances r y - have the 
form 

lim M'k, .k, [Vij. Rfc) -> (27r)" 3 exp(z k„- ■ r y - + i K fc • R fe ) 

-Rfc— >oo 3 



exp r±ia i3 - ln(fcy r y - ± k y - ■ r y -) j , V (r y -, 
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The '+' and '— ' signs refer to outgoing and incoming boundary conditions, respec- 
tively. Similarly to the two-body problem, the Sommerfeld-parameter are given 
by ' 

n„ = (12) 

The asymptotic state (|II]) is a straightforward generalization of the two-body asymp- 
totic given by Eq.(|2|) to three-body systems. However, unlike the situation in two- 
body scattering, in three-body systems other types of asymptotic are present where, 
in a certain set (r^, R&), one Jacobi coordinate tends to infinity whereas the other 
coordinate remains finite ||. The asymptotic states ( ]TT| ) serve as the boundary 
conditions that have to be satisfied by the scattering solutions of the Schrodinger 
equation. The derivation of these solutions is a delicate task and will be the subject 
of the remainder of this section. 

The general approach here is to consider the three-body system as the subsume 
of three non- interacting two-body subsystems ||. Since we know the appropriate 
coordinates for each of these two-body subsystems (the parabolic coordinates) we 
formulate the three-body problem in a similar coordinate frame with 

{£fc = r ij ± • r ij}> e ijk 0; j > i,k e [1, 3], (13) 



where ky denote the directions of the momenta k^. Since we are dealing with a six- 
dimensional problem three other independent coordinates are needed in addition 
to (0). To make a reasonable choice for these remaining coordinates we remark 
that, usually, the momenta are determined experimentally, i. e. they can be 
considered as the laboratory-fixed coordinates. In fact it can be shown that the 
coordinates (|T3|) are related to the Euler angles. Thus, it is advantageous to choose 
body-fixed coordinates. Those are conveniently chosen as 

{6 = nj}, e ijk ^ 0; j>i,ke [4, 6]. (14) 



Upon a mathematical analysis it can be shown that the coordinates (0,0) are lin- 
early independent || except for some singular points where the Jacobi determinant 
vanishes. The main task is now to rewrite the three-body Hamiltonian in the coor- 
dinates (0,|TJ). To this end it is useful to factor out the trivial plane- wave part [as 
done in Eq.@] by making the ansatz 

v i'k k, Rfe) = N exp(i r i:j ■ + % R fc • K fc ) v|/ k Ki (ry, Rfe)- (15) 
Inserting the ansatz (15) into the Schrodinger Eq.fllOD leads to the equation 

3 



— A Yi] + — A Rk + 2i (—kij ■ V r .. + —K fc • V Rfe ) - 2 — 
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#(r o -,R fc ) = 0. 



(16) 



In terms of the coordinates (p~3| , [T^) Eq. (|16|) casts 

//^k .<, ((:• • • • = [ H par + H in + tf mix ]^ 3 , Kfc (6, • • -,6s) = . (17) 

The operator H par is differential in the parabolic coordinates £1,2,3 only whereas Hi nt 
acts on the internal degrees of freedom £4,5,6- The mixing term H m i X arises from the 
off-diagonal elements of the metric tensor and plays the role of a rotational coupling 
in a hyperspherical treatment. 

The essential point is that the differential operators H par and H int are exactly sep- 
arable in the coordinates £i... 3 and £ 4 .. . 6 , respectively, for they can be written as || 
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H par = Y,H^] [H £j ,Ht i ] = Q;Vi,je {1,2,3}, (18) 

i=i 

and 
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H int =J2 H ti' \- H tv H ^ = 0; V 2, j G {4, 5, 6}, (19) 
where 

= [dfy £j + iki m £j d^. — fii m ; 

ejim^O, 3 e {1,2,3}. (20) 
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The operator H mix = H — H par — H int derives from the expression 
H mix ■= EU=i K V ^) • + (V Rfc £ u ) • (Vr^,)} 



(21) 



(22) 



(23) 



(24) 



Noting that H^ ,j = 1,2,3 is simply the Schrodinger operator for the two-body 
scattering rewritten in parabolic coordinates (after factoring out the plane-wave 
part), one arrives immediately, as a consequence of Eq.(^), at an expression for the 
three-body wave function as a product of three two-body continuum waves with the 
correct boundary conditions ([ll]) a ^ large inter-particle separations. This result is 
valid if the contributions of H int and H mix are negligible as compared to H par , which 
is in fact the case for large interparticle separations |J or at high particles' energies. 

It should be noted that this same result can deduced in a Jacobi coordinate 
system however the operators H par , H int and H mix have a much more complex 
representation in the Jacobi coordinates (cf. Ref. 0). 

With decomposing the total Hamiltonian in H par , Hi nt and H m i X we achieved a 
result similar to that obtained by Pines [I] and co-worker for the interacting electron 
gas: The system is decomposed into a long-range and a short range components 
described respectively by H par and H int . As clear from Eq. flTTD , the eigenfunctions 
of H par have an oscillatory asymptotic behaviour whereas the eigenstates of Hi nt 
decay for large interparticle distances |L4]| . The mixing term H m i X couples the short- 
range to the long-range modes of the system. 

The analytical structure of the Eqs.flH 



deserves several remarks: 



The total potential is contained in the operator H par , as can be seen from 
Eqs.flSHp. Thus, the eigenstates of H par treats the total potential in an exact 
manner. This means on the other hand that the operators H int and H mix are 
parts of the kinetic energy operator. This situation is to be contrasted with 
other treatments |i~7H25|l of the three-body problem in regions of the space 
space where the potential is smooth, e.g. near a saddle point. In this case 
one usually expands the potential around the fix point and accounts for the 
kinetic energy in an exact manner. 

In Eq.(^) the total potential appears as a sum of three two-body potentials. 
It should be stressed, that this splitting is arbitrary, since the dynamics is 
controlled by the total potential. I.e., any other splitting that leaves the total 
potential invariant is equally justified. This fact we will use below for the 
construction of three-body states. For large inter-particle separation the op- 
erators Hi n t and H mix are negligible as compared to H par and the splitting 



of the total potential as done in Eqs. (pUP becomes unique. This means, for 
large particles' separation the three-body dynamics is controlled by sequential 
two-body scattering events. 



The momentum vectors kjj enter the Schrodinger equation via the asymptotic 
boundary conditions. Thus, their physical meaning, as two-body relative mo- 
menta, is restricted to the asymptotic region of large inter-particle distances. 
The consequence of this conclusion is that, in general, any combination or 
functional form of the momenta kjj is legitimate as long as the total energy 
is conserved and the boundary conditions are fulfilled (the energies and the 
wave vectors are liked via a parabolic dispersion relation). This fact has been 
employed in Ref. || to constructed three-body wave functions with position- 



dependent momenta kj ? - and in Ref. [26| to account for off-shell transitions. 



The separability of the operators (|2l| - |23|) may be used to deduce representa- 



tions of three-body states jnj that diagonalize simultaneously H par and H int . 
It should be noted however, that generally the operator H mix , which has to be 
neglected in this case, falls off with distance as fast as H int . 

As well-known, each separability of a system implies a related conserved quan- 
tity. In the present case we can only speak of an approximate separability and 
hence of approximate conserved quantum numbers. 

If we discard H int and H mix in favor of H par , which is justified for 7^ 
0, k G [1, 3] (i.e. for a large or for a high two particle momentum kij), the 
three-body good quantum numbers are related to those in a two-body system 
in parabolic coordinates. The latter are the two-body energy, the eigenvalue of 
the component of the Lenz-Runge operator along a quantization axis z and the 
eigenvalue of the component of the angular momentum operator along z. In 
our case the quantization axis z is given by the linear momentum direction k^. 

In Ref. || the three-body problem has been formulated in hyperspherical- 
parabolic coordinates. In this case the operator Hi nt takes on the form of the 
grand angular momentum operator. This observation is useful to expose the 
relevant angular momentum quantum numbers in case H m i X can be neglected. 

In Ref. H 27|] the three-body system has been expressed in the coordinates 
r/j = £f, j = 1,2,3 and fjj = j = 1,2,3. This is the direct exten- 
sion of the parabolic coordinates for the body problem (cf. Section 3) to the 
three-body problem. From a physical point of view this choice is not quite 
suited, for scattering states are sufficiently quantified by outgoing or incom- 
ing wave boundary conditions (in contrast to standing waves, such as bound 
states whose representation requires a combination of incoming and outgoing 
waves). Therefore, to account for the boundary conditions in scattering prob- 
lems, either the coordinates rjj or fjj are needed. The appropriate choice of the 
remaining three coordinates should be made on the basis of the form of the 
forces governing the three-body system. In the present case where external 
fields are absent we have chosen = r^, k = 4, 5, 6 as the natural coordinates 
adopted to the potential energy operator. 



3.1. Coupling the short and the long-range dynamics 



In the preceding sections we pointed out that the eigenstates of H par can be de- 
duced analytically. These eigenf unctions, even though are well defined in the entire 
Hilbert space, constitute a justifiable approximation to the exact three-body state 
in the asymptotic region only (e.g. in the region of larger inter-particle separation 
or for higher energies). This fact is important when it comes to evaluating reaction 
amplitudes, for such amplitudes involve the many-body scattering state in the entire 
Hilbert space. Therefore, an adequate description of the short-range dynamics may 
be necessary, in particular in cases where the contributions to the matrix elements 
of the transition amplitudes originates from the internal region, i.e. when the re- 
action takes place at small interparticle distances. Nevertheless the eigenstates of 
H par can, and have been used for the calculations of transition matrix elements, as 
for example done below. In this case the justification of this doing must go beyond 
the asymptotic correctness argument. 

In this section we seek three-body wave functions that diagonalize, in addition 
to H par , parts of H int and H mix . 

One method that turned out to be particularly effective for this purpose relies 
on the observations: a) In a three-body system the form of the two-body potentials 
Zij/rij are generally irrelevant, as long as the total potential is conserved, c) To keep 
the mathematical structure of the operators (|T8| , |2~0|) unchanged and to introduce a 
splitting of the total potential while maintaining the total potential's rotational 
invariance one can assume the strength of the individual two-body interactions, 
characterized by Zij, to be dependent on £4,5,6- This means we introduce position 
dependent product charges as 



z ij(£,4, £5, Ce) ; (25) 



with 
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To obtain the many- body potentials := Zij/r^ we express them as a linear mixing 
of the isolated two-body interactions Vij := z^/r^, i. e. 



(27) 



where A(£a, £5,^) is a 3 x 3 matrix. The matrix elements are then determined 
according to 1) the properties of the total potential surface, 2) to reproduce the 
correct asymptotic of the three-body states and 3) in a way that minimizes Hint and 
H mix . It should be stressed that the procedure until this stage is exact. It is merely 
a splitting of the total potential that leaves this potential and hence the three-body 
Schrodinger equation unchanged. 



3.2. An electron pair in the field of a positive ion 

To be specific let us demonstrate the method for the case of two electrons moving 
in the Coulomb field of a residual ions. This brings about some simplifications 



since the ion can be considered infinitely heavy as compared to the electron mass. 
Traditionally, the electrons are labeled a and h and their positions and momenta 
with respect to the residual ion are respectively called r a , r b and k a , k b . Adopting 
this notation, the eigenstate of the operator H par reads: 

^k a ,k b (£l- 6 ) = 1^1 073a, 1, -iK fl ) 

1F1 (i(3 b , 1, -ik b £ 2 ) 
iFi (i/3 ab , 1, -ik ab £ 3 ) . 

(28) 

Here we denoted the relative electron-electron momentum by k ab = |(k a — k b ) 
whereas i-Fi[a, b, x] stands for the confluent hypergeometric function and /3j, j G 
{a, b, ab} are the Sommerfeld parameters 

P 1 = ^l, je{a,b,ab}, (29) 

with Vj being the velocities corresponding to the momenta k, and Zj, j G {a, b, ab} 
are the electrons-ion and electron-electron effective product charges, respectively. 
The form of zj is still to be determined. 

Below the functions Zj^.-s) are given that preserve the total potential, possess 
the correct three-body asymptotic and incorporate features of the many-particle 
motion at the complete fragmentation threshold, namely along the saddle point of 
the total potential, the so-called Wannier ridge flTl — p5| . Since Zj are assumed to 
depend on the internal coordinates only (£4, £5,^6) the wave function fl2"Bp is still an 
eigenstate of the long-range Hamiltonian H par (given by Eq.|20|) . In physical terms 
it can be said that the effect of the short-range part of the Hamiltonian is to modify 
dynamically the coupling strength of the isolated two particle system (%•). 

To ensure the invariance of the Schrodinger equation under the introduction of 
the product charges Zj{^...^) the three conditions 



E- 



:,(C;M,) - Z + — + —,je{a,b,ab}, (30) 



Tj r a r b r ab 

have to be fulfilled (r ab is the electron-electron relative coordinate and z is the charge 
of the residual ion). The wave functions containing Zj must be compatible with the 
three-body asymptotic boundary conditions. These are specified by the shape and 
size of the triangle formed by the three particles (two electrons and the ion): I. e., 
the derived wave function must be, to a leading order, an asymptotic solution of the 
three-body Schrodinger equation when the aforementioned triangle tends to a line 
(two particles are close to each other and far away from the third particle) or in the 
case where, for an arbitrary shape, the size of this triangle becomes infinite. The 
latter limit implies that all interparticle coordinates r a ^^ ab must grow with the same 
order, otherwise we eventually fall back to the limit of the three-particle triangle 
being reduced to a line 0, as described above. In addition we require the Wannier 
threshold law for double electron escape to be reproduced when the derived wave 
functions are used for the evaluation of the matrix elements. 

All of the above conditions are sufficient to determine zj and thus the wave func- 



tion (28). This wave function is called "dynamically screened three-body Coulomb 



wave function" ^dssc- This is because this wave function consists formally of three 



Coulomb waves where the short-range dynamics enters as a dynamical screening of 
the strength of the two-body interaction. 

The applicability of the wave function \1/ds3C to scattering reactions is hampered 
by the involved functional dependence leading to complications in the numerical 
determination of the normalization and of the scattering matrix elements. Further- 
more, the incorporation of the three-body scattering dynamics at shorter distances 
brings about intrinsic practical disadvantages as compared to an approach where 
Zj are constant (this approach is usually called three-body Coulomb wave method 
, 3C). Namely, the construction of ^dssc nas to be individually undertaken for 
given charge and mass states of the specific three particle system at hand. This is 
comprehensible since properties of the total potential are inherent to the particular 
three-body system under investigation. 



As shown in Ref. |28| the normalization of the 3C wave function is readily de- 
termined from the asymptotic flux. This procedure has not been accessible in the 
case of \1/ds3C due to the position dependence of Zj. 

To overcome this difficulty (and that associated with the six-dimensional nu- 
merical integration when evaluating transition matrix elements) we note that the 
position dependence of Zj(r a ,r b ,r ab ) occurs (due to dimensionality considerations) 
through ratios of the interparticle distances. Thus, this dependence can be converted 
into velocity dependence by assuming that 

^oc^. (31) 

r 3 V 3 

The proportionality constant in Eq. (|3l|) could be of an arbitrary functional depen- 
dence. It should be emphasized that the approximation fl3"T| ) is not a classical one, 
i. e. it is not assumed that the particles' motions proceed along classical trajectories 
[conversely, if the motion were classically free, Eq. fl3"T| ) holds]. It merely means 
that the total potential is exactly diagonalized in the phase space where Eq. (|3~ID is 
satisfied, as readily deduced from Eq. (|30|) . 

Eq. (pip renders possible the normalization of \I/ds3C since in this case we obtain 



z 
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Zj(k a} kb,k a b) and the arguments used in Ref. [£8|] can be repeated to deduce 



for the normalization N the expression 
N = I[ N 3 J e {a,b,ba} 

3 

Nj = exp[-P j (k a ,k b ,k ba )Ti/2]T[l-ip j (k a ,k b ,k ba )]. (32) 



Here T(x) is the Gamma function. The velocity-dependent product charges [29| 
have the form 



Zba(y a ,Vb) = [1 - (/ gf a bl ] a b * (33) 

^a(Va,V 6 ) = -1 + (1-Z ba 



(34) 

v 1+a 

^b(v a , v 6 ) = -1 + (1 - z ba ) b — . 

K + V b) V ab 

(35) 



The functions occurring in Eqs. ( |33| , |34D are defined as (v , v fe are the electrons' 
velocities and v ab = v a — v&) 

3 + cos 2 4a v a 

j := , tana = — (do) 

J v b 

(37) 
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6i := 38 



2vaV b cos(6 ab /2) 

g 2 (-0.5 + p) (39) 
E (40) 



£ + 0.5' 

where is being measured in atomic units and p is the Wannier index (the value 
of p depends on the residual ion charge value, the numerical value of p for a unity 
charge of the residual ion is p = 1.127). The interelectronic relative angle 6 ab 
is given by 9 ab := cos -1 v a • v^. In case of higher excess energies (E ^> 1) it is 
readily verified that a — > 1 [Eq. fl4"0|) 1 and all modifications of the charges ([33"l - [3~5l ) 



which are due to incorporating the Wannier threshold law become irrelevant. The 



charges ( p3| - p5| ) reduce then to those given in Ref. || with Eq. (|3l| ) being applied. 
From the functional forms of the charges (|33}-P^|) it is clear that when two particles 
approach each other (in velocity space) they experience their full two-body Coulomb 
interactions, whereas the third one 'sees' a net charge equal to the sum of the charges 
of the two close particles. 



3.3. Applications to atomic scattering problems 

In this section we assess the analytical methods developed above by performing 
a numerical evaluation of many-body scattering amplitudes. The reaction we are 
considering here is the electron-impact ionization of atomic hydrogen. In the final 
channel of this collision process two interacting electrons move in the double con- 
tinuum of a residual ion. Hence a correlated three-body wave function is needed to 
represent this state. For this wave function we employ the approximate expressions 
given in the preceding sections. The initial state consists of an incoming single- 
particle wave that represents the projectile electron and a bound state of atomic 
hydrogen. 

The complete information on this reaction is obtained by measuring the coinci- 
dence rate for the emission of two continuum electrons with specified wave vectors, 
i.e. the energies E a , E b and the emission solid angles Q a , Q b of the two electrons are 
determined for a given incident energy of the projectile electron. Due to energy con- 
servation it suffices to determine the energy of one of the electrons. Therefore, one 
measures in this way a triply differential cross section (TDCS), i.e. a cross section 
differential Q a , Q b and E b . 

If the spin of the electrons is not resolved, the TDCS is a statistically weighted 
average of singlet and triplet scattering cross sections 

TDCS(k a , k b ) = c Q|T S | 2 + -IT*] 2 ) (41) 

where k« is the momentum of the incident projectile and c = (27r) 4 (k a k b )/ki. The 
singlet T s and triplet transition matrix elements T* derive from the corresponding 



transition operators T s and T*, where 



T s = (l + V ab )T fi (k a ,k b ) 

T = (l-V ab )T fi (KM)- (42) 

The action of the exchange operator V ab on the operator 7}j is given by T ab Tfi(k a , k 6 ) - 
Tfi(k b , k a ). The prior representation of 7}j(k a , k&) is given by 

T /f (k a ,k 6 ) = (*|F i |$ k4 ) ■ (43) 
The wave function \1/ is obtained from Eq. ( p8|) as 

*k a ,k a = iVexpi(k a ■ r a + k 6 • r 6 )^ kajka . 

The three-body system in the initial channel is described by l^iq)- Assuming |$kj) 
to be the asymptotic initial-state, i. e. (r , r&l^&i^) is a product of an incoming 
plane wave representing the incident projectile electron and an undistorted ls-state 
of atomic hydrogen, the perturbation operator V\ occurring in Eq. (fi3| ) is given by 
l/|r„ — r&| — l/r a (which is the part of the total Hamiltonian not diagonalized by 
|$ki))- In what follows we choose the the incident direction kj. The final 

state electrons are detected in a coplanar geometry, i. e. kj • (k a x k b ) = 0. The z 
axis lies along the direction perpendicular to the scattering plane, i. e. parallel to 
k a x kft. The polar and azimuthal angles of the vector k a (k&) are denoted by 6 a , <ft a 
(6 b ,4> b ), respectively. In the coplanar geometry considered here the polar angles 
are fixed to 9 a — 7r/2 = Q b . In the calculation of the DS3C model we employ the 
approximation (|3lD and use the product charges ( |3~3"l - [35l ) . If we use the unit matrix 



for the transformation (E7I), i.e. if we assume A = I, the three-body wave function 
reduces the eigenfunction of the asymptotic part H par of the Hamiltonian without 
any coupling to the internal region. This wave function is commonly known as 
the 3C wave function |3],|28|]. In addition we compare the results of the analytical 



methods presented here with those of the convergent close coupling method (CCC). 
This is a purely numerical method that attempts at evaluating exactly the transition 
matrix elements fully numerically. In Fig.l the angular distribution of one of the 
electrons is shown for two fixed angular positions of the other electron. The two 
electrons are ejected with equal energies E a = E b = 6.8 eV. As clear from Fig. 1 the 
effect of the coupling to the internal region is very important, since the results of 
the 3C model that neglects the short-range dynamics are at clear variance with the 
experiment. The differences between the CCC method and the experiments are still 
the subject of current research. The main advantage of analytical methods is that 
they allow an insight into the origin of the structures observed in the cross sections. 
An extensive analysis underlying this statement has been carried out in Ref.( ||31|| ) 
where the main peaks in Fig. 1 have been assigned to certain sequence of collisions 
between the participating particles. 

At higher incident energies the discrepancies between the DS3C and the 3C re- 
sults disappear and both of those models (as well as the CCC calculations) are in 



overall agreement with the experiments |36|. From this situation one can conclude 



that at higher energies the short-range parts of the Hamiltonian (H int and H mix ) 
are of less importance, for they have been neglected in the 3C model whereas the 
DS3C theory accounts for them via the dynamical screening (we note the asymptotic 
region is reached for large k^k, i- e - f° r large momenta the distance does not need 
to be very large). 




Figure 1: The fully differential cross section for the electron- impact ionization of 
atomic hydrogen in the co-planar, symmetric energy-sharing geometry. 
The incident energy is = 27.2 eV. One ejected electron is detected at 
a fixed angle <3? a with respect to the incident direction [ $ a = 30° (upper 
panel) and <3> a = 45° (lower pannel)]. The angular distribution of the 
other emitted electron is measured. The emission angle of this electron 
with respect to the incident direction is denoted by Both electrons have 
the same energy, namely E a = E/, = 6.8 eV. Full squares are experimental 
data Ref. |3D|,^T|]. The solid thick lines show the predictions of the DS3C 
theory employing the matrix A (cf. Eq.(^)) whereas the dotted curves 
indicate the results of the 3C theory, i.e. when using A = 1. Thin solid 
lines are the full numerical calculations using the convergent close coupling 
method (CCC). 
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Figure 2: The single differential cross section for the electron-impact ionization of 
atomic hydrogen as function of the ratio E a /E where E a is the energy of 
one of the final-state electrons and E = E a + E b is the total excess energy 
which is chosen as E = 200 meV. The use of the \1/ds3C approximation 
yields the solid curve whereas the dashed curve represents the results when 
employing the 3C model (A — I). The 3C results have been multiplied by 
a factor of 10 13 for a better shape comparison. 



This means in physical terms that the two electrons attains their asymptotic mo- 
menta swiftly without much of scattering from intermediate states whose behaviour 
is determined mainly by H int and H mix . 

Integrating over all emission angles of the two electrons we end up with a single 
differential cross section depending on the energy of one of the electrons. Since the 
energy of the other electron is then determined via the energy-conservation law, 
the single differential cross section has to be symmetric with respect to the point 
where both electrons have the same energy. Fig. 2 shows the results for the single 
differential cross sections as calculated within the DS3C method along with the 
calculations within the 3C method. The excess energy is very low (200 meV). For 
small excess energies the Wannier theory, which relies on phase space arguments, 
predicts a flat energy distribution between the electrons, i.e. a flat single-differential 
cross section. This prediction has been substantiated by full numerical calculations 
[57 . As seen in Fig. 2 the DS3C predicts a flat energy sharing between the electrons 



close to the complete fragmentation threshold, in contrast to the 3C results which 
are strongly peaked around the equal energy- sharing configuration. This deviation 
of the 3C results from those of the Wannier theory is not surprising since in the 
Wannier approach one expands the potential around a saddle point (accounting for 
terms up to a fourth order) and neglects higher order terms while the kinetic energy 
is treated fully. In contrast the 3C model neglects the short-range part of the kinetic 
energy. Obviously it is this part which is most important for the Wannier mode and 
the resulting predictions. 

Sampling over the energy sharing between the two electrons, i.e. integrating the 
single differential cross section shown in Fig. 2, one obtains the total cross section 
as function of the excess energy E = E a + E\, (or equivalently as function of the 
incident energy Ei). Close to the three-body break-up threshold the total cross 
section cr(E) for two continuum electrons receding from a charged ion has been 
investigated by Wannier |[7J using a classical analysis. Wannier [T7J pointed out 



that the excess-energy functional dependence of the total ionization cross section at 
the three-particle fragmentation threshold can be deduced from the volume of the 
phase space available for double escape of the two electrons. For the present case 
of atomic hydrogen Wannier deduced the threshold law cr(E) oc E 1A2r . Since then 



an immense amount of theoretical and experimental studies ( e.g. ||18|--p5|, p8H41||) 
using quite different approaches have been carried out which basically confirm the 
Wannier-threshold law. 

The Wannier treatment predicts the scaling behaviour of the cross section cr(E), but 
it does not provide any information about the magnitude of cr(E). That the magni- 
tude is a very sensitive quantity is illustrated by the behaviour of the cross section 
in the independent Coulomb particle model which is obtained in our case by switch- 
ing off the interaction between the two electrons in the final channel. In this case 
the cross section reveals a linear dependence on the excess energy, cr(E) oc E [0. 
Although the latter dependence of cr(E) does not deviate much from the Wannier 
threshold law {cr(E) oc E 1A2r ) the absolute value of cr(E) within the independent 
Coulomb particle model is largely overestimated [compare Fig. 3]. If we employ the 
wave function \I/ds3C; with the dynamical product charges described in the preceding 
sections we end up with results in good accord with the experimental measurements 
(cf. Fig. 3). The absolute magnitude of the total cross section is satisfactorily repro- 
duced when the DS3C model is employed. To examine the analytical behaviour of 
(j {Ei) calculated using \Pds3C we plot in Fig. 3(b) the quantity a(E)/E 1A27 . Accord- 




Figure 3: The total cross section cr(E) for the electron-impact ionization of atomic 
hydrogen as a function of the excess energy E. The solid (long dashed) 
curve shows the results for u{E) when treating the two continuum electrons 
according to the DS3C theory (3C model) whereas the dashed dotted 
curve denotes the results of the independent Coulomb particles model (see 
text). Results of the CCC method are also included (short dashed curve). 
Experimental data are due to Shah et al. | 3~2|] . The inset in the panel 
(a) shows the results of the 3C theory on a logarithmic scale. In the 
upper panel (b) the quantity cr(E)/ E 1A27 is depicted as a function of E as 
evaluated using \I/ds3C- 
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Figure 4: The spin asymmetry [Eq.(|44])l in the total ionization cross section u{E) for 
the electron-impact ionization of atomic hydrogen. Results of the hidden 



crossing theory |25| (long dashed curve) and those of the CCC method |53" | 
(solid curve) are depicted. Treating the final state according to the DS3C 
theory (3C theory) yields for the spin asymmetry the results shown by 
the short-dashed (dot-dashed) curve. The experimental data are due to 
Fletcher et al. |34|] (full squares) and Crowe et al. [|35]| (open circles). 



ing to the Wannier-threshold law the latter quantity should be a constant function 
of E and gives the absolute value of the cross section. As seen in Fig. 3(b) the Wan- 
nier threshold law is in fact reproduced by the cross section results of the \I/ds3C 
within a range of E 6 [0, 0.5e^/]. For E > 0.5 eV the analytical dependence of <j{E) 
evaluated with \I/ds3C slowly deviates from the Wannier threshold law. When using 
the 3C method for the description of the two escaping electrons (A = I in Eq. (|27|) ) 
we obtain an analytical behaviour for a(E) which is not compatible with the Wan- 
nier theory. The absolute value for the total cross section is as well not reproduced 
by the 3C model, for the reasons discussed above. Also included in Fig. 3 are the 
results of the convergent-close coupling method, CCC, fl3~3"| . The results of the CCC 



are in good agreement with the experimental <j(E) for higher energies |p5| , however, 
close to threshold the evaluation of <j(E) is limited by the computational resources 
as an ever increasing number of pseudo states is needed to achieve convergence. 
In addition to the magnitude of the cross section, the spin asymmetry, A, offers a 
further way of probing the dynamical properties of the electron-impact ionization of 
atomic systems. The spin asymmetry A is defined as 

„, v a s (E)-a t (E) , . 

where a s and <r* are the total ionization cross sections for singlet and triplet scatter- 
ing, respectively. The Wannier theory for threshold ionization predicts a constant 
value of A with increasing excess energy but provides no information on the numer- 



ical value of A j|3 . Measurements of A at threshold reveals a slightly positive slope 



of the spin asymmetry with increasing excess energies ||44]|. In Fig. 4 the results for 



A are shown in the case where the two-electron continuum final state is described 
by the 3C theory and by \I/ds30 Also depicted in Fig. 4 are the results of the CCC 



approach []33] and the method using hidden-crossing theory |25| . Although all theo- 
ries, except for the 3C model, are in reasonable agreement with experimental finding 
the positive slope of A at threshold is not reproduced. 

Neglect of the short-ranged part of the Hamiltonian H in and H mix , i.e. using the 
3C model, results in a completely wrong behaviour of the calculated spin asymmetry. 
With increasing excess energy the inner region of the Hilbert space becomes of less 
importance for the present reaction and the results of the 3C method become more 
and more in better agreement with the experimental data. 

We note here that since the spin asymmetry is a ratio of cross sections it is 
expected that the spin asymmetry is rather sensitive to the detailed of the radial part 
of the wave functions. From the agreement between the experiment and the DS3C 
theory observed in Fig. 4 we conclude that the radial part of \I/ds3C is well behaved 
at lower excess energies and that the short-range parts of the total Hamiltonian 
Hi n t and H m i X plays a dominate role at lower energies, as far as the value and the 
behaviour of the spin asymmetry are concerned. 



4. Correlated states of N charged particles 



In the preceding sections we considered the two and three-body continuum spectrum. 
Unfortunately, the curvilinear coordinate system (|l3"||[4Tj used for the three-body 
problem does not have a straightforward generalization to the N body case. There- 
fore we will treat the problem of N charged particles at energies above the complete 
fragmentation threshold within a reference frame spanned by a set of Jacobi coordi- 
nates. The problem becomes more transparent if we consider N — 1 particles with 
equal masses (and with charges Zj, j £ [1, N — 1]) that move in the Coulomb field 
of a residual massive charge z. The mass M of the charge z is assumed to be much 
larger than m (M ^> m). Neglecting terms of the order m/M the centre-of-mass 
system and the laboratory frame of reference can be chosen to be identical. The 
non-relativistic time-independent Schrodinger equation of the iV-body system can 
then be formulated in the relative-coordinate representation as 



N N 
j=l '3 i4 1 13 



E 



',3 

j>i=l 



*(ri, 



r N ) 



(45) 



where rj is the position of particle j with respect to the residual charge z 



-1-3 •' 



— Tj denotes the relative coordinate between particles i and j 



and 
The kinetic 



energy operator H has the form (in the limit m/M — > 0) H = — ^ £=1 A£/2m 
where is the Laplacian with respect to the coordinate r^. We note here that 
for a system of general masses the problem is complicated by an additional mass- 
polarization term which arises in Eq. (^5j). Upon introduction of iV-body Jacobi 
coordinates, Hq becomes diagonal, however, the potential terms acquire a much 
more complex form. 

Assuming the continuum particles to escape with relative asymptotic momenta k, 
(with respect to the charge z) it has been suggested in Ref. 0, due to unpub- 
lished work by Redmond, that for large interparticle distances the wave function 



^(rx, • • - , r N ) takes on the form 

N N 

lim ¥( rij ... jr „) - (2ny m / 2 HUrs)Mrs) H^jfaj) > ( 46 ) 



i,j = l 

j>i 



where the functions £j( r j)> i[>j(rj), i>ij{ T ij) are defined as 

0( r j) : = ex P( ik i- r i) (47) 
■= exp ^piaij ln(fcj r,- ± k, • r,) , (48) 



V>«(r«) := exp 



Tiaij ln(% ry ± ky ■ r^) . (49) 



The + and — signs refer to outgoing and incoming boundary conditions, respectively, 
and ky is the momentum conjugate to r^, i. e. k^ := (kj — kf)/2. The Sommerfeld- 
parameters a,, c% are given by 

oca = — -, a.j = — -. (50) 

In Eq. flSUP i>j denotes the velocity of particle j relative to the residual charge 
whereas Vy := Vi — vj. In this work we restrict the considerations to outgoing- wave 
boundary conditions. The treatment of incoming-wave boundary conditions runs 
along the same lines. The total energy of the system E is given by 



N , 2 

E = Y, E u where E l = -i-. (51) 



2m 

To derive asymptotic scattering states in the limit of large inter-particle separations 
and their propagations to finite distances we assume for ^(i - !, • • • , r N ) the ansatz 

#(r l5 • • • , v N ) = jV$/(ri, • ■ • , ^^//(ri, • ■ ■ , t n )x(ti, ■ ■ ■ ,r N ) (52) 

where $7-7- are appropriately chosen functions, AT is a normalization constant 
and x( r i> • • • , r^v) is a function of an arbitrary form. The function is chosen to 
describe the motion of iV-independent Coulomb particles moving in the field of the 
charge z at the total energy E, i. e. is determined by the differential equation 

N \ 

H + J2 — -E)<S>i(r 1 ,---,rN) = 0. (53) 
3=1 rj / 

Since we are interested in scattering solutions with outgoing-wave boundary condi- 
tions which describe iV-particles escaping with asymptotic momenta k,, j e [1, AT], 
it is appropriate to factor out the plane-wave part and write for 

N 

S/fa,...,!^) =¥/(r 1 ,---,r JV ) JJ&fa). ( 54 ) 

3=1 



Upon substitution of the ansatz (m) into Eq. Q53|) it is readily concluded that Eq. 
(|53|) is completely separable and the regular solution can be written in closed 
form 

N 

$i(r 1 ,---,r N ) = Y[Z j (r j )tp j (T j ) (55) 
3=1 



where (fjivj) is a confluent-hypergeometric function in the notation of Ref. |46| 



= lFl[<*j, 1, -OkfTj + kj ■ Tj)]. 



(56) 



The function <&j describes the motion of the continuum particles in the extreme case 
of very strong coupling to the residual ion, i. e. \zzj\ ^> \zjZi\] Vi,j G [l,iV]. In order 
to incorporate the other extreme case of strong correlations among the continuum 
particles (\zjZi\ ^> \zzj\\ V i,j G [1, N]) we choose $// to possess the form 



$ /7 (r 1; • • • , r N ) = $//(r 1; , r N ) £(i 
with 

AT 

$//(ri,---,rAr) := | [ ^(r^) 



(57) 



(58) 



j>i=i 



where <Pij(rij) '■= l-Fifa^-, 1, —i(kijrij + k^ ■ r^)]. It is straightforward to show that 
the expression ^-(r^-) YliLi Ci( r i) solves for the Schrodinger Eq. (|45|) in the case 
of extreme correlations between particle i and particle j , i. e. \zzi\ <C \ziZj\ 3> 



V l,m,n 7^ i, j. In terms of differential equations this means 

A 7 



#0 + 



Z%Zj 



'J 



e ) r^i^ru 

i=i 



0. 



(59) 



It should be stressed, however, that the function fl5"?D does not solve for Eq. ([45|) 
in case of weak coupling to the residual ion (z —>■ 0), but otherwise comparable 
strength of correlations between the continuum particles. This is due to the fact 
that two-body subsystems formed by the continuum particles are coupled to each 
other. To derive an expression for this coupling term we note first that 



m— 1 



N 



N 



N 



A m $/j = A m y ;m 17 (fjj + ^ A m (p mn Yl <Pij + An, me[l,JV](60) 



i=i 



n=m+l 



3>i 



where the differential operator A m has the form 



A, 



m—l 

1=1 

m—l 

+£ 

i=i 



N 



(V ' m^l m ) 



N 

+ E 

n=m+l 



( ^ ^ V ' m^Prnn) 
n=m+l 

m—l 

( ^ ^ ^ ' m^Psm) 
l^s=l 

N 



N 



n ^ 



j>i 

N 

n ^ 



(^ ' mSPmri) 



t=m+l 



N 



J J <ftj m G [l,iV]. 



(61) 



To obtain the differential operator which couples the two-body subsystems in ab- 
sence of the charge z we neglect in (^5|) the interactions between the residual charge 
and the continuum particles (z = 0) and substitute the function ( |57"1) into Eq. 



Making use of the relation fl6"0|) it is straightforward, however cumbersome, to show 
that the coupling term which prevent separability has the form 



iv 



A=J2 A m- (62) 



m=l 



Eqs. ( |6T| , |62"D warrant commented upon: The term A m is a mixing operator. It 
couples an individual two-body subsystem formed by two continuum particles to 
all other two-body subsystems formed by the continuum particles in absence of the 
residual ion. Hence it is clear that all the terms in the sum (|6l|) vanishes for the 
case of three-body system since in this case only one two-body system does exist in 
the field of the residual charge. The second remark concerns the structure of A m 
and hence A. From Eq. flB"0| ) it is evident that the remainder term ( |51~D is part of the 
kinetic energy operator. Thus it is expected that, under certain circumstances, this 
term has a finite range which indicates that asymptotic separability, in the sense 
specified below, does exist for many-body continuum Coulomb systems. In fact, as 
the functional form of (pijfaj) is known the term A can be calculated explicitely 
which will be done below. 

Now with and $// have been determined, the exact wave function ( |52"1) is given 




by the expression x( r i; ' ' ' j v n)- Upon substitution of the expressions (|57|j55|) into 
the ansatz (|52|) and inserting in the Schrodinger equation (f45|) a differential equation 
for the determination of x( r i; ' ' ' j t n) is derived 

N 

+ (V* In $ 7 ) ■ (V, In $ /7 ) ] + X (ri, • • - , v N ) = . (63) 

From the derivation of the functions and [Eqs. (|53l , |57|) 1 it is clear that all long- 
range two-body Coulomb interactions have been already diagonalized by and 
because the total potential is exactly treated by these wave functions. Hence, the 
function x, to be determined here, contains information on many-body couplings, 
which are, under certain conditions (see below), of finite range. To explicitely show 
that, and due to flux arguments we write the function \ in the form 

N 

X (n, , r N ) = H f( rj )[l - /(r x , • • • , r N )] (64) 
i=i 

where /(r 1; • • • , rj\r) is a function of an arbitrary structure. Inserting the form (|64T) 
into Eq. (|63D we arrive, after much differential analysis, at the inhomogeneous dif- 
ferential equation 

#o-X^(ln$/ + ky$ n )+iki]-Vi\f + K{l-f) = (65) 
t=i J 

where the inhomogeneous term 1Z is given by 

N ( m-1 N 

TZ:=Y^ \ (V m ln$/) ■ (V m ln$//) + ^ ( V m In <p lm ) ■ (V m In <p mp ) 

m=l \ 1=1 p=m+l 



^ lit ±111, _L 

+ - / l(Vw m Plm) ■ (V OT In 

^ am ) 

1=1 s^l 

N N ^ 

) ■ ( V m In 

n=m+l n^q=m-\-l ) 



1=1 s^l 

N ~\ 

It is the inhomogeneous term 7Z which contains the coupling between all individ- 
ual two-particle subsystems. For example the first term in Eq. flBTf ) describes the 
coupling of a two-body subsystems formed by particles i and j to all two-body sub- 
systems formed by the individual continuum particles and the residual ion. The 
second term originates from (|62l) and, as explained above, is a measure for the cou- 
pling among two-body subsystems of the continuum particles (in absence of z). To 
these couplings to be negligible the norm of the term 1Z must be small. For get some 
insight into the functional form of 7Z, given by fl66|), we note that 

V e \n$ I = a e k e F e {r e ) (67) 

where 

„ f v [1 + iae, 2, -i{k e r e + k e - r e ) ] f 

1F1 [ia e , 1, -H^ r e + k e - r t ) \ 

In addition we remark that 

N m-1 

V m ln$// = ^ V m ln V m In V9 im 

n=m+l Z=l 

iV m-1 

^ ^ ^ rank mn^ ran (j^mn) ^ ^ ^Im^lm^ lrn(? Im) (69) 



n=m+l i=l 



where 



F / r \ ._ 1^1 [1 + jag, 2, -ijkjj + kjj ■ Tjj) ] f 

:/•', [iay, 1, -/(/.;, r« + k. • r,,) ] ^ + ' 1 ' 

Thus, the behavior of the coupling term 7Z is controlled by the generalized functions 
Fy(ry), Fj|(r/) since Eq. (|66|) can be written in the form 



N 



TZ ■ ^ S < ^m^mF m (r m i 



m=l 



N m-1 



^ ^ O^mnkmnP mn (j^mn) ^ ^ ^sm^smFsin(r s 



n=m+l s=l 



m—1 N ^ m—1 m—1 



^ ^ ^ ^ (-^Im^mpklmkmpP lm ' F mp ^ ^ ^ ^ ^ImP^ svnHlmH srnf ' F sm 

Z=l p=m+l Z=l s^Z 

1 N N \ 

n=m+l n^q=m+l J 

The simplest approximation is to neglect the term 1Z altogether. In this case the 
function / = solves for the equation (|65|) . Then, the solution of Eq. ( |45| ) takes on 
the approximate form 

N 

V(r 1 ,---,T N )t*N Y[ Awnfam). (72) 

m>l,j=l 



Thus, the justification of the approximation ([72]) reduces to the validity of neglecting 
the inhomogeneous term ([7T|) . One region in which this term can be disregarded 
is the asymptotic region of large inter-particle separations. This is immediately 
deduced from the asymptotic behavior of the generalized functions Fj^r^), F;(r;) 
which dictate the asymptotic properties of TZ, as readily concluded from Eq. (|7I|). 



From the asymptotic expansion of the hypergeometric functions we infer that 



lim \F ij (r ij/ 



ij 



kjj • (kj 



r ij) r ij 



O (\kij + kjj ■ Yij\ 2 ) . 



(73) 



A asymptotic relation similar to Eq. ([73]) holds for F^(r;). It should be noted that 
the functions Fj^r^), F;(r/) have to be considered in a distributive (operator) sense 



which means that, asymptotically, only terms of Fy , F; which fall off faster than the 
Coulomb potentials can be disregarded. Since TZ is essentially a sum of products of 
Fij, Fi the expression TZ is of finite range, in the sense that it diminishes faster than 
the Coulomb potential in the asymptotic regime, only in the case where all particles 
are far apart from each other, i. e. 



lim TZ 



O 



n- 



\ki ri + k r ri\ 2 ) 



V j>i,le[l,N}. (74) 



Therefore, in the limit (0), the term TZ can be asymptotically neglected and the 
approximation (|72|) is justified. In fact, it is straightforward to show that the wave 
function ([T^) tends to the asymptotic form fl4"6| ) in the limit of large inter-particle 
separations which proves the assumption made in Ref. 0. However, if two particles 
are close together, regardless of whether all other particles are well separated, the 
coupling term is of infinite range, as seen from Eqs. (|73|j71|) . In this case the relation 
(|74l) does not hold. Consequently, the wave function (|72|) is not an exact asymptotic 
eigenfunction of the total Hamiltonian in this limit. It is important to note that 
the limit Eq. fl74l) is energy dependent. With increasing momenta of the escaping 
particles the asymptotic region, i. e. the limit Eq. (|71|), is reached faster. In other 
words, at a certain inter-particle separations, the remainder term TZ, which has 
been neglected to arrive at the approximate form (|T2D, diminishes with increasing 
velocities of the emerging particles. In this sense the approximation leading to the 
wave function (|7^) is a high energy approximation. 



4.1. The two-body cusp conditions 



In the preceding section it has been shown that the approximation (|72|) is, to leading 
order, exact for large particles separation. In addition, it is concluded below that 
this function exhibits a behavior compatible with equation fl45|) at all two-body 
coalescence points ry — > 0,r L — > 0, j > i,l G [l,iV]. To guarantee regular behavior 
of the wave function at these collision points, at which the corresponding Coulomb 
two-body potential is divergent, the solution \I/(ri, • • ■ , rjv) of Eq. (EH) must satisfy 



the Kato cusp conditions |45, 47[] (provided the solution does not vanish at these 
points). At a collision point r% — > these conditions are 



d^(n,---,r N ) 

d r, ; 



0. 



V {n/ rj ) -> 0, {n/r lm ) -»■ 0; m > I, i ± j 6 [1, N] (75) 



The quantity ^(rj., • • • , r N ) is the wave function ^(r 1 , ■ ■ ■ , r N ) averaged over a sphere 
of small radius < 1 around the singularity rj = 0. A relation similar to Eq. (f75|) 
holds in the case of the coalescence points — > 0. To prove that the wave function 
(|72|) does satisfy the conditions (|75D we linearize \?(ri, • • • , r N ) around = and 
average over a sphere of small radius 1 to arrive at 



N 



^(ri, • • • , r N ) = M D{vi) Y\ £j¥>j(rj)pim(rim), e«m 



(76) 
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where 
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l + «j hi n . 



r\ dcos9 



1 + ikiCosO + ttj/cj rAl + cost 



To arrive at Eq. (|77D one takes the z— axes as kj and define cos# 
Eqs. (|77| , |7"6"1) it is obvious that 



(77) 
From 
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In deriving Eq. ([73) we made use of the fact that in the limit (r^/r^ — > 0) the 
distance tends to r 3 -. The proof that the wave function ([7^) fulfills the cusp 



conditions at the collision points of two continuum particles 



0) runs along 



the same lines. Finally, we remark that the wave function (|72|) is not compatible 
with the expansion of the exact solution of the Schrodinger equation ( f45[) at the 
three-body collision points (e.g. r$ — > and rj — > 0,j ^ i) since in this case the 
exact wave function is known to satisfy a Fock expansion |15| in the coordinate 



p := y (r? + r|) which contains, in addition to powers in p, logarithmic terms in p 

) possesses a regular power-series expansion around 



whereas the wave function 
0. 



Ti — > and rj 



4.2. Normalization 



The knowledge of the normalization factor M of the wave function (|72|) is imperative 
for the evaluation of scattering amplitudes using the wave function ( |T2| ) as a rep- 
resentation of scattering states. In principle, M is derived from a 3iV-dimensional 
integral over the norm of the function (|72|) which, for large N, is an inaccessible 
task. Thus for the determination of TV we resort to the requirement that the flux 
through an asymptotic manifold defined by a constant large inter-particle separa- 
tions should be the same in the case of the wave function (|72|) and a normalized 
plane-wave representation of the scattering state, i. e. 



Jpw — J* 

where the plane-wave flux is given by 



(79) 
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In Eq. ( j30|) the total gradient V := Yli=i nas been introduced. To evaluate the 
flux generated by the wave function ([72]) we note that, by taking advantage of Eqs. 



\p9), we can write for the total gradient of the wave function 
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where F mn is given by F mn {p mn . The decisive point now is that since we are consider- 
ing the flux at large interparticle distances only the first term of Eq. (^1|) is relevant. 
This is readily deduced from Eqs. (|68|j70|) which state that all other terms in Eq. 
flS~T|), except for the first term, can be neglected asymptotically. Note in this context 
that terms in the wave function which are asymptotically of the order 0(l/rj, l/rj m ) 
correspond to parts of the Hamiltonian falling off faster than the Coulomb poten- 
tials and hence can be disregarded in the asymptotic regime. Now making use of 
the asymptotic expansion of the confluent hypergeometric function [I5| and taking 
leading order in the interparticle distances the flux can de deduced 
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where T(x) is the Gamma function. From Eqs. ([79| , |80|J82] ) it follows that 
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For two charged particles moving in the field of a heavy nucleus the wave function 
(|72|) with the normalization, given by Eq. (|83|), simplifies to the three-body wave 
function that has been discussed in the previous section. 



5. Green function theory of finite correlated systems 

In the preceding sections we investigated the two, three and iV-body correlated 
scattering states. With increasing number of particles the treatment becomes more 
complex and a methodology different from the wave function technique is more 
appropriate. A method which is widely used in theoretical physics is the Green 
function approach which we will follow up in this section. 

For a canonical ensemble, we seek a non-perturbative method which allows to 
distribute systematically the total energy between the potential and the kinetic 
energy parts. This is achieved by the development of an incremental method in 
which the N correlated particle system is mapped exactly onto a set of systems in 
which only N — M particles are interacting (M e [1, N — 2]), i.e. in which the 



potential energy part is damped. (In contrast to re-normalization group theory 
we do not reduce the strength of interactions, but the number of them). This is 
particularly interesting from a thermodynamic point of view since for a number of 
thermodynamic properties the kinetic energy contributions can be separated out 
from the potential energy parts, as shown in the next section for the internal energy. 
By virtue of the present method the potential energy part is systematically reduced. 

For a formal development let us consider a nonrelativistic system consisting of 
N interacting particles. We assume the total potential to be of the class = 
Ylj>i=i v ij without any further specification of the individual potentials Vij. For 
three-body potentials the development of the theory proceeds along the same lines. 



The potential satisfies the recurrence relations 



N 



um = F^5>r'' (84) 

£ <" 2 '. 3 * k > (85) 
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Uj ~ N-3 
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where vJf* ^ is the total potential of a system of N — 1 interacting particles in which 



3 



the j particle is missing, i.e. in terms of the physical pair potentials v mn , one can 
write uf 1 ' = Em>n=i v mn, m^j^n. 

The fundamental quantity that describes the microscopic properties of the N 
body quantum system is the Green operator G^ N ' which is the resolvent of the total 
Hamiltonian. It can be deduced from the Lippmann Schwinger equation G^' = 
Co + GoU^G^ where Go is the Green operator of the non interacting N body 
system. An equivalent approach to determine the dynamical behavior of a system is 
to derive the respective transition operator which satisfies the integral equation 
T {N) = jj(N) + u(N)G T^ N \ These integral equations for G (JV) and provide a 
natural framework for perturbative treatments. However, for N > 3 the application 
of the above Lippmann Schwinger equations (and those for the state vectors) is 
hampered by mainly two difficulties: 1.) as shown in Refs. |48|,|49j the Lippmann 



Schwinger equations for the state vectors do not have a unique solution, and 2.) 
as shown by Faddeev [p>0|-|52"|] the kernel of these integral equations K = GqU^ 
is not a square integrable operator for iV > 3, i.e. the norm \\K\\ = [Tr(KK^)] 1 / 2 
is not square integrable. The kernel K is also not compact. The reason for this 
drawback is the occurrence of the so-called disconnected diagrams where one of the 
N particles is a spectator, i.e. not correlated with the other N — 1 particles. For the 



three-body problem Faddeev ||50| , |5T| suggested alternative integral equations with 
square integrable kernel. 

Our aim here is twofold: (a) We would like to derive non-perturbative integral 
equations that treat all iV particles on equal footing and are free from disconnected 
diagrams, (b) These equations should allow to obtain, in a computationally acces- 
sible manner, the solution of the correlated N body problem from the solution when 
only N — M particles are interacting (where M G [1,N — 2]). 

According to the decomposition (0), the integral equation for the transition 
operator can be written as 
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T ( N) = y T (N-i) (g6) 
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Here we introduced the scaled potential 
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The transition operator of the system, when N — 1 particles are interacting via 



the scaled potential u- 



is 



With this relation Eq . (p7|) can be reformulated as 
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Eq.(B3) can be expressed in a matrix form as follows 
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The kernel [K^ ^] is a matrix operator and is given by 
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(90) 



From Eq. (|85D it is clear that tj N ^ can also be expressed in terms of the transition 
operators of the system where only N — 2 particles are interacting: 

N-l 

f (N-l) _ T (N-2) 
k+3 



The operators T^ N ^ are deduced from Eq.([S9"l) with N being replaced by N — 1. 



From the relation = G + G T^G o we conclude that the Green operator 
of the interacting A" particle system has the form 
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The operators G^ ^ are related to the Green operators l ' 



which only A" 
given via 
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(92) 



\ N / 



where [K^" 1 )] = G [K^ N ^]Go \ From Eqs.(|8|j9|) we conclude that if the 
Green operator of the interacting N — 1 body system is known (from other analyt- 
ical or numerical procedures, e.g. from an effective field method, such as density 
functional theory) the Green operator of the N particles can then be deduced by 
solving a set of N linear, coupled integral equations (namely Eqs.( P^|J92] )). Accord- 
ing to the above equations, if only the solution of the N — M problem is known 
where M G [1, AT — 2] we have to perform a hierarchy of calculations starting by 
obtaining the solution for the N — M + 1 problem and repeating the procedure to 
reach the solution of the A^ body problem. 

At first sight the kernels of Eqs. (|89| , |92"D appear to have disconnected diagrams 
since they contain transition operators of systems where only A^ — 1 particles are 
interacting and one particle is free (disconnected). It is, however, straightforward 
to show that any iteration of these kernels is free of disconnected terms (the discon- 
nected terms occurs only in the off-diagonal elements of [K^ - ^ and [K^ - ^). For 
N = 3 the present scheme reduces to the well-established Faddeev equations. As for 
the functional structure of the Eqs. ( |55| , P^ ) we remark that for the solution of the 
A^ particle problem we need the (off-shell) transition operators of the A^ — 1 subsys- 
tem. The interaction potentials do not appear in this formulation (in contrast to the 
Lippmann Schwinger approach). On the other hand the (on-shell) transition matrix 
elements can be determined experimentally. This fact becomes valuable when the 
potentials are not known. 



5.1. Application to four-body systems 

Over the years a substantial body of knowledge on the three-particle problem has 
been accumulated. In contrast, theoretical studies on the four-body problem are 
still scare due to computational limitations whereas an impressive amount of exper- 
imental data is already available [SB-pT/j]. Thus, it is desirable to apply the above 



procedure to the four-body system and to express its solution in terms of known 
solutions of the three-body problem. For A^ = 4 the first iteration of Eq. (^3) yields 

4 

G w = X>f - 3G . (93) 



Here is the Green operator of the system where only three particles are inter- 
acting and can be taken from other numerical or analytical studies. This means, to 
a first order, methods treating the correlated three-body problem can be extended 
to deal with the four-body case using Eq.(|93[). We note that for the case of non- 
interacting system reduces to = G and hence Eq.flDlf) reduces to = Gq, 
as expected. 

The Green function encompasses the complete spectrum of the many-body sys- 
tem, i.e. the wave function approach can be retrieved from the Green function. For 
example, Eq. (^) leads to an expression for the four-body state vector in the form 

l* (4) > = l^) + I^S) + l^i) + l^)-3|0i 4 e ) e ). (94) 

Here \ifj^l) is the state vector of the system in which the three particles i, j and k 
are interacting whereas |0^ e ) is the state vector of the non- interacting four-body 
system. The state vectors \ipif k ) can be approximated by Eq.(|28|) or by the other 
procedures discussed in the preceding section on the three-body problem. 

Since the state vector ( |94"|) is expressed as a sum of correlated three-body states, 
the evaluation of the four-body transition matrix elements for a specific reaction 
simplifies considerably. In addition, the spectral properties of a many-body inter- 
acting system can be obtained in a straightforward way from those for systems with 
a reduced number of interactions, for in this case the matrix elements of the total 
Green functions are expressed as sums of matrix elements of reduced Green func- 
tions, as evident from Eq. (|93|) . This spectral feature can be exploited to study the 
thermodynamical properties of finite correlated systems. 



6. Thermodynamics and phase transitions of 
interacting finite systems 

To investigate the thermodynamical properties of N interacting particle system we 
remark that at the critical point divergent thermodynamical quantities, such as the 
specific heat Cy are obtained as a derivative with respect to the inverse temperature 
(3 of the logarithm of the canonical partition function Z(f3), 

C v =p 2 dl\nZ(P) = f((3,Z((3))/Z({3). 

Here / is some analytical function and for the Boltzmann constant we assume k — 1. 

Therefore divergences in the thermodynamic quantities, which signify phase tran- 
sitions are connected to the zero points of Z(/3). These zero points are generally 
complex valued. Therefore an analytical continuation of Z(f3) to complex tempera- 
tures is needed. 

The connection between the phase transitions and the complex zero points of the 
grand canonical partition function have been uncovered by Yang and Lee |5£| . In this 
case one seeks an analytical continuation of the fugacity z = exp(/3/z) (here \x is the 
chemical potential) to the complex plane z — > §l(z) +iQ(z). In the thermodynamic 
limit the zero points condense to lines. The transition points are the crossing points 
of these lines with the real fugacity axis. 

Grossman et al. generalized the concept of Yang and Lee to the canonical 
ensemble. In this case the inverse temperature (3 = 1/T is continued analytically to 
(3 = 3t{f3) + i^s(f3). The phase transitions are then the crossings of the zero points 



line of Z(0) with the real /3 axis. The advantage here is that a classification of the 
phase transitions can be given in terms of how the zero-points line do cross the real 
3 axis |EU . 



The crucial point is that in the thermodynamic limit iV — > oo, V — > oo and 
v = V/N < oo (V is the volume) the zero points approach, to an infinitesimal 
small distance the real axis. For this reason, the characteristic phase-transition 
divergences appear in the thermodynamical quantities. For finite systems Z(f3) has 
only finite zero points that do not approach infinitely close the real axis. Therefore, 
the thermodynamic quantities show smooth peaks rather then divergences. The 
positions and widths of these peaks can be obtained from the real and imaginary 
parts of the zero points laying closest to the real axis |K) . 



To apply this method to correlated finite systems we need a representation of the 
canonical partition function that can then be continued analytically to the complex 
temperature plane. 

The canonical partition function of a correlated system can be expressed in terms 
of the many-body Green function as 

Z(f3) = J dE n(E) e~ m . (95) 

Here Q(E) is the density of states which is related to imaginary part of the trace of 
G^via 

ME) = --QTrGW(E). (96) 
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From the Green function expansion Eq. Q9~ip we deduce 
Z(J3) = J dETrG {N \E) e~ pE 
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To a first order Zj = Zj N ^ — Z where Z^" 1 ^ is the partition function of a system 
in which only N — 1 particles are interacting. For the applications of the Grossmann 
method let us remark that Z(j3) is an integral function and can be expressed in a 
polynomial form. Recalling the analytical properties of meromorphic functions one 
can write Z((3) in terms of its complex zero points as 

Z'(0) / 3\ p 

z(3) = z{oy^w TT e^. (100) 



7. Conclusions 



In conclusions we reviewed recent advances in the analytical treatment of corre- 
lated few charged-particle systems. Starting from the two-body (Kepler) problem 
we discussed the mathematical and physical properties of a three-body system above 
the complete fragmentation threshold. Analytical approximate expressions for the 
three-body wave functions have been derived and employed for the calculations of 
transition matrix elements in atomic scattering experiments. The discussion has 
been extended to the case of N charged particles in the continuum where N par- 
ticle wave functions have been derived and their mathematical features have been 
exposed. For the description of the complete spectrum of a general finite system 
we discussed a Green function method that maps the N interacting particle system 
onto a system with a less number of interactions. A brief account has been given 
on how this method can be used to investigate thermodynamic properties of finite 
systems. The application of the Green function method for the calculations of scat- 
tering amplitudes have been presented in Ref. The Green function method 



presented here have been extended recently to deal with the scattering of correlated 
systems from ordered and disordered potentials Q . 
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